/A/ "3 "2^ «* ^3SC~>0 



NASA Technical Paper 3552 


Closed-Form Evaluation of Mutual Coupling in 
a Planar Array of Circular Apertures 


M. C. Bailey 


April 1996 




NASA Technical Paper 3552 



Closed-Form Evaluation of Mutual Coupling in 
a Planar Array of Circular Apertures 


M. C. Bailey 

Langley Research Center • Hampton, Virginia 


National Aeronautics and Space Administration 
Langley Research Center • Hampton, Virginia 23681-0001 


April 1996 



Available electronically at the following URL address: 

Printed copies available from the following: 

NASA Center for AeroSpace Information 
800 Elkridge Landing Road 
Linthicum Heights, MD 21090-2934 
(301) 621-0390 


http://techreports.larc.nasa.gov/ltrs/ltrs.html 


National Technical Information Service (NTIS) 
5285 Port Royal Road 
Springfield, VA 22161-2171 
(703) 487-4650 



Abstract 


The integral expression for the mutual admittance between circu- 
lar apertures in a planar array is evaluated in closed form. Very good 
accuracy is realized when compared with values that were obtained by 
numerical integration. Utilization of this closed-form expression, for 
all element pairs that are separated by more than one element spacing, 
yields extremely accurate results and significantly reduces the compu- 
tation time that is required to analyze the performance of a large elec- 
tronically scanning antenna array. 


Introduction 

The wide flexibility available in the design of antenna arrays is very useful in applications where 
factors such as beam shaping, side-lobe control, and rapid beam steering are of prime consideration; 
however, the implementation of a good design can become quite complicated as a result of the effects of 
mutual interaction between closely spaced radiating elements. These interactions are evident as (1) a 
distortion of the radiation pattern, (2) an element-driving impedance that varies as the array is phased to 
point the beam in different directions, and (3) a polarization variation with scan angle in an array with 
elements that can support more than one sense of polarization. The degree to which the interelement 
coupling affects the performance of the array will depend upon the element type, the polarization and 
excitation of each element, the geometry of the array, and the surrounding environment. To accurately 
model the effects of mutual interelement coupling in the design of a phased-array antenna, the analysis 
must include all these factors. 

Since interelement coupling in phased arrays is a near-field phenomenon, an accurate analytical 
model is generally formulated such that the resulting expression involves either a single or a double 
integration in the spectral domain. This integral formulation can readily be evaluated numerically with 
the aid of high-speed computers; however, the computation time can still become prohibitively large. 
This substantially increased computation time for large arrays is primarily a result of the need to calcu- 
late the mutual coupling between all possible pair combinations of the array and the highly oscillatory 
nature of the integrand to be evaluated, which oscillates more rapidly and converges more slowly as the 
separation between element pairs increases. 

The focus of this paper is to illustrate a technique for developing an accurate closed-form evaluation 
of the integral formulation for mutual coupling between circular-aperture elements in a planar array. In 
particular, the final results in this paper are limited to identical circular elements whose aperture fields 
are restricted to that of the dominant mode of a circular waveguide of the same cross section; however, 
the approach is applicable to other aperture fields that can be represented in Bessel-function form. 
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quantities defined in equation (34) 
electromagnetic wave frequency, rad/sec 


Analysis 

The general analytical formulation for the interelement mutual coupling in planar arrays of arbitrary 
apertures has been developed. (See ref. 1 .) The effects of mutual coupling are determined by computing 
the self and mutual admittances among all the elements of the array to form a complex admittance 
matrix. This admittance matrix is then operated on to determine the complex scattering matrix for the 
array. The scattering matrix represents the relationship between the amplitudes and phases of all the 
feed-waveguide modal fields that are incident on and reflected from the radiating apertures. This com- 
plex scattering matrix allows one to completely characterize the performance of the antenna for any 
amplitude and phase excitation. 

The mutual admittance between the ith and yth apertures of the array can be determined from the 
reaction between the electric field of the ith aperture and the magnetic field of the jth aperture. In 
general, 


Y ■■ = 
y V. V 


:J/[ ] 

' s, 


E (i \x, y,0)xH U \x,y,0) 


z dx dy 


(1) 


where E*-^(x,y,0) is the vector electric field in the ith aperture with all others long circuited, and 
H^(x,y,0) is the vector magnetic field that would exist at the ith aperture with all apertures short 
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circuited except the jth. The fields are evaluated at the aperture plane (z = 0). However, for purposes of 
this analysis (which will become obvious), equation (1) is rewritten as 


Y ij = ^JJ[ e(, W 0) >< hm H °W,z) 


■ z dx dy 


( 2 ) 


In the Fourier spectral domain, equation (2) can be expressed in an equivalent form (ref. 1) as 


y </ = 7S 
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z dk r dk ,, 

x y 


(3) 


By performing the vector multiplications and by indicating that the integrations are to be performed 
before taking the limit, the admittance expression is rewritten as 
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1 J 
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(4) 


In the spectral domain, the transverse components of the transformed electric and magnetic fields are 
related to the transformed solutions to the wave equations as 
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(5) 


where the primes on A and F denote differentiation with respect to z. The admittance expression can 
then be rewritten as 
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( 6 ) 


When the first term is multiplied and divided by A\-k x ,—ky, 0), and the second term is multiplied and 
divided by F(-k x ,-ky,0), equation (6) becomes 
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Solving equation (5) independently for the zth and jth apertures and substituting into equation (7) 
yields the following general expression for the mutual admittance between two apertures whose 
electric-field distributions are known: 


,2 fc o 
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( 8 ) 


where a change of variables has been made such that k x = k 0 fi cos a and k y - /c 0 p sin a. The solutions to 
the free-space wave equations in the Fourier transform domain yield 


rfe 0 e r (0)A(a, p,zri 
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-jk 0 zj 1 -P 2 j 

_ jA\ a, P, 0) . 
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For circular apertures, whose field distributions are those of the circular waveguide modes, the integra- 
tion on a in equation (8) can be readily evaluated in terms of Bessel functions. The mutual admittance 
expression for circular apertures then reduces to a single integration on (3. 

To demonstrate the method of evaluating the closed-form expression for mutual admittance, the 
remainder of this report is limited to identical circular apertures with only the domi n ant transverse elec- 
tric mode (TEn) aperture fields. The mathematical development for other circular waveguide modal 
fields and unequal size apertures would proceed in the same manner. Assuming identical TEn mode 
circular apertures, equation (8) can be written as 
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( 10 ) 


Although the two apertures are identical in size and excitation, they may be polarized differently with 
respect to each other (as denoted by the relative polarization angle <j> p ). The geometry for the two circu- 
lar apertures is illustrated in figure 1. 

To evaluate equation (10) in closed form, the semiconvergent series of Hankel (ref. 2, pp. 137 
and 138) is first utilized to express the Bessel functions (with arguments k 0 $R) in a series form as 

exp H) 

= v — :( r v (~^oP^) ex P(^oP /? i~^ T v (^oP i? i ex P(~^oP /? )i (H) 
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where 


00 f . n n 1 

T v (k 0 $R) = 1 + x { ~ j) - „ n t 4v2 - ( 2m - ^ 
n = l inKM o m n m = l 


( 12 ) 


and where v is the order of the Bessel function. (For the special case of TEn mode, v is either 0 or 2.) 
When equation (11) is substituted into equation (10) and P is substituted for — P in terms containing 
T v (-k 0 $R), the integration on P in equation (10) can be extended over the limi ts of (-°° to °°) as 


where the following additional quantities have been introduced for convenience: 
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Expressing z and R in spherical coordinates (z = r cos 0 and /? = rsin0), equation (13) can be 
rewritten as 


12 


\ liT C^expfj^Um { ~7== P g(P)exp[;>/(p)] 

K ^0 H k 0 r V 4; e _ >5 [7iin0 j -~ 


( 21 ) 


where 



( 22 ) 
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The integral in equation (21) is now in a form that can be readily evaluated, for large values of r, by 
the saddle-point method (ref. 3, pp. 305 to 307) as 


g(P)exp[jr/(P)] dp = g(P 0 ) ex P[J r /(Po)5 f~~7 

» — oo k —11 
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The primes on /and g denote differentiation with respect to P, and p 0 is the saddle point as determined 
from 


/'(P 0 ) = 0 

/ (Pq)^® 


(24) 


Therefore, 


P 0 = sin 9 


(25) 


The modification in equation (2) allowed the saddle point to be determined as defined in equation (24); 
thus, the integral evaluation can be performed. 

The evaluation of the integral in equation (23) requires taking partial derivatives (up to the fourth 
order) with respect to the integration variable P and evaluating these derivatives at the saddle point Po- 
The quantities of interest are 
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and 
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The primes on W\ and W 2 denote differentiation with respect to p. Taking derivatives of g with respect 
to p, setting P = p 0> and dropping the functional notations yields 
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All primes in equations (30) and (31) denote derivatives with respect to p evaluated at P = p 0 . Substitut- 
ing equations (26) to (31) into equation (23) and evaluating equation (21) at 0 = tc/ 2 yields (after consid- 
erable algebraic manipulation) 
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Retaining only the terms in HR, HR 2 , and HR 2 , the mutual-admittance expression becomes 
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Equation (35) can be used to calculate the complex mutual admittance between a pair of identical circu- 
lar apertures excited in the dominant TE U mode; however, the self admittance of a single aperture must 
still be calculated by numerical integration (ref. 1). 

The self admittance can be obtained from the integral form of the mutual admittance. 
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(36) 


by making the two apertures coincident (i,e., by setting both the polarization angle and the aperture 
spacing to zero). The aperture self admittance is then given by 
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(37) 


In a large array, the self admittance for each element and the mutual admittance for each pair com- 
bination is calculated; and the complex scattering matrix for the array can then be determined from the 
complex normalized admittance matrix by the following matrix relationship: 


with 
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rs,i 


■fV 


r 5 i/] + t% i 
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(i = j ) 1 
V*j) J 



(38) 


(39) 


(40) 


where 5 ( y is the complex coupling coefficient from aperture j to aperture i, F ( is the characteristic admit- 
tance for the waveguide mode that excites aperture i, and [ ] _1 denotes matrix inversion. 

The scattering matrix, in conjunction with the array excitation coefficients, contains the neces- 
sary information for describing the performance of a scanning phased-array antenna, including all the 


10 



interactions between array elements. The scattering matrix gives the relationship between the incident 
and reflected waveguide modal fields for all the elements of the array as follows: 


b \ = Sjjaj +S\2 a 2 + ^13 a 3 + "• + *^1 N a N 
b 2 = ^21 a l + ^ 22 a 2 + ^23 a 3 + + ^2N a N 

b 3 = ^31 a l + ^32 a 2 + ^33 a 3 + "* + ^3N a N f 


b N ~ S N l a l + S N2 a 2 + S N3 a 3 + ”* + S NN a N) 


The equivalent voltage and current for the ith aperture is then expressed as 

v i = a . + b i | 
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Therefore, the active admittance of the ith element in the array is 
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and the active reflection coefficient of the ith element is given by 
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Assuming a constant incident power source, the radiated beam from the planar array can be scanned 
to the angle (9 0 , <t>o) by producing a progressive phase shift across the array as follows: 


a ( = a t exp[->A: o (A: ; cos(|) o + y ! .sin<l) o )sin0 o ] 


(45) 


where x t and y, are the Cartesian coordinates of the center of the ith aperture. The variation with beam 
scan for the reflection coefficient of a particular element in the array can be calculated by substituting 
equation (45) into equation (44) with the scattering coefficients obtained from equation (38). 


Results 

The closed-form expression (eq. (35)) for mutual admittance was validated by comparing results 
with those obtained from the numerical integration of equation (36). The scattering-matrix results were 
used for comparison; in all cases, the self admittance was obtained by numerical integration of equa- 
tion (37). All computations were performed on an MS-DOS 80486DX2-66 desktop computer using the 
code CWG (ref. 4); this code was modified to include an option for the closed-form mutual-admittance 
evaluation. 

Two sets of data were used for verification. The first data set was the mutual coupling between two 
apertures, and the second set was the active reflection coefficient for a large scanning array whose grid 
geometry is illustrated in figure 2. 

The first set of data was obtained for the purpose of establishing the range of applicability by com- 
paring the mutual coupling between two identical apertures for various sizes, orientations, and spacings. 
These results (plotted in figs. 3 and 4) indicate that very accurate results are consistently obtained for 
center-to-center spacings greater than the classic far-field distance of 2 D 2 /Xq; in many cases, accurate 
results are obtainable when the apertures are almost touching. 
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The second set of data was calculated for the active reflection coefficient of the center element in a 
scanning array. The mutual admittance between all element pairs for a 721-element array (array diame- 
ter of 20 wavelengths) was calculated by numerical integration and by the closed-form expression. The 
calculated reflection coefficient for the center element is plotted against beam scan angle in the two 
principal planes in figures 5 to 8. The results of these calculations, using numerical integration for all 
element pairs, were used as the basis for evaluation. When the closed-form expression was used for all 
element pairs, a small discrepancy was observed, as shown in figures 5 and 6. However, when the 
admittance matrix was modified to use the numerical integration values for the closest neighbor pairs 
(i.e., elements spaced 0.714 Xq) and the closed-form expression for all others, the results were extremely 
accurate, as shown in figures 7 and 8. As a result of the uniform grid geometry of the array, numerical 
integration of only two mutual-ad mi ttance values were required in order to obtain extremely accurate 
results. Also, as noted in figures 7 and 8, utilization of the closed-form expression resulted in a consid- 
erable reduction in computational time required to fill the admittance matrix. 

Conclusion 

An approach for obtaining a closed-form expression for the mutual admittance between elements in 
a planar array of apertures was presented. The closed-form expression for circular apertures was devel- 
oped and compared with results obtained by numerical integration. The judicious use of the closed-form 
expression, in conjunction with the integral form of the mutual admittance, provides an antenna design 
and analysis tool that produces extremely accurate results with a significant reduction in computational 
time for large phased arrays. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
December 11, 1995 
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Figure 1. Geometry of dominant transverse electric mode (TE n ) excited circular apertures. 
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Figure 2. Geometry for equilateral triangular grid array. 
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Figure 5. Reflection coefficient versus H-plane scan angle for center element of 721-element array. 
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Figure 6. Reflection coefficient versus E-plane scan angle for center element of 721-element array. 




17 




Reflection coefficient amplitude Phase, deg 




Beam scan angle, deg 

Figure 7. Reflection coefficient versus H-plane scan angle for center element of 721 -element array. 
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Figure 8. Reflection coefficient versus E-plane scan angle for center element of 721 -element array. 
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